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Abstract 

By a completely independent analytic calculation made in general covariant 
gauge we confirm the known result for the massless 0(a^) correction to Rhad 
and check explicitly its gauge independence. We extend the calculations to 
include contributions due to a hypothetical existence of a colour octet of (light) 
neutral fermions known also under the name of gluinos. In numerical form and 
for five active flavours we found (in MS scheme) that R(s) — y[l + "s/ 71 " + 
(1.409 - 0.346n 5 )(a s /vr) 2 + (-12.805 - 3.006n g - 0.0466n?) (a s /vr) 3 ] where n g 
is the number of such colour octets. 
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1 Introduction 



The simplest inclusive reaction involving quarks — their production through a decay of a 
heavy virtual photon or a Z boson — is a process of fundamental importance for QCD as 
the theory of strong interactions. Firstly, it provides us with a beautiful confirmation of 
the existence of colour Secondly, the precision measurements of F(Z — > hadrons) have 
developed into an important experimental tool for a reliable determination of a s (Mz) @. 
Thirdly, by comparing a s [Mz) with the value of a s as obtained from the measurements 
of R{s) at lower energies one could in principle directly test the evolution of the strong 
coupling constant ||. 

Impressive wealth of theoretical results about R{s) is now available in the QCD frame- 
work (see e.g. a recent review M for a detailed discussion). In the massless approximation, 
valid in the high energy limit, corrections to R = cr(e + e~ — ► hadrons)/cr(e + e~ — ► 
have been calculated up to order |5], |j . For precision measurements the dominant mass 
corrections can be included through an expansion in m 2 /s. Terms up to order alm 2 /s 
and a 2 s my s 2 || are available at present, providing an acceptable approximation from the 
high energy region down to intermediate energy values. At order a 2 even full account of 
the quark mass effects has been recently made on the basis of a semi-analytic approach 
||. Among massless calculations the one in the order otl is probably most involved. Let 
us discuss it separately. 

The results of the reevaluation of massless next-next-to-leading correction to R{s) 
(first computed erroneously in Ref. [|10| because of bugs in a computer program) were 
published in Refs. pi Q five years ago. Both calculations have produced identical results 
and are much alike in many respects: 

• the same theoretical tools as well as computer programs have been used; 

• the simplest gauge condition - the Feynman one - has been employed; 

• An important subclass of all diagrams (which includes, in fact, the most complicated 
ones) - the QED type diagrams (that is those not containing three and four-gluon 
vertices) - had been computed first and published in Ref. JTT]] by a collaboration of 
authors of Refs. || |J. 

A golden rule well-known among multi-loop people says that a result of a multi-loop 
calculation can be trusted and considered as the result only if it is confirmed by an 
independent calculation preferably made by a different group and with the use of the 
general covariant gauge. Therefore, in view of the importance of the 0(a^) correction for 
both theory and phenomenology it is necessary to check those by a really independent 
calculation. 

In the present paper we report about such an attempt. Our work confirms the re- 
sults of Refs. 0, |(J and explicitly demonstrate their gauge independence. We have also 
computed extra diagrams which include virtual loops of Majorana fermions transforming 
as an octet with respect to the colour SU C (3) group. Note that the case of Majorana 
fermions considered as gluinos, that is as superpartners for gluons ||12j| , has some phe- 



nomenological relevance at present. On one side, light gluino are not completely excluded 



by the current experimental data [13[ [14], 15] . On the other side, overall consistency of 
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various determinations of a s (Mz) seems to improve if the light gluino with mass of a few 
GeV does exist |T5], [T7|, |I§. 



2 Preliminaries 

From purely theoretical point of view R{s) is an extremely suitable object to deal with 
within pQCD. This is because all relevant information is contained in the current corre- 
lation function 

n^(g) = J dxe^(0\T[ j™(x)j™(0) ]|0) = (-g, u q 2 + q^M-q 2 ) , (1) 

with the hadronic EM current = Y^fQfi^fl^f, and Qf being the EM charge of the 
quark /. The optical theorem relates the inclusive cross-section and thus the function 
R(s) to the discontinuity of IT in the complex plane 

R(s) = —1mU{-8-iS). (2) 

4:71 

Experimental e + e~ data are taken in the physical regime of timelike momentum transfer 
q 2 > 0. This region is influenced by threshold and bound state effects which make the 
use of perturbative QCD questionable. However, perturbative QCD is strictly applicable 
for large spacelike momenta (q 2 = —Q 2 < 0), since this region is far away from non- 



perturbative effects due to hadron thresholds and resonance effects |LJ| . Therefore, reliable 
theoretical predictions can be made for II (Q 2 ) with Q 2 > 0. To compare theoretical 
predictions and experimental results for time-like momenta, one has to perform suitable 



averaging procedures [|0|, |2]J . For large positive s one may appeal to the experimentally 
observed smoothness of R as a function of s and to the absence of any conceivable non- 
perturbative contribution. 

The renormalization mode of the polarization operator IT(Q 2 ) reads (see, e.g. Ref. M) 

n(Q 2 //i 2 ,a s ) = Z™ + II (Q 2 ,a°),, (3) 

where a s = g 2 /(47r) is the strong coupling constant. Within the MS scheme |22) (here 
and below we are using a convenient combination a s = a s /7i) 



z cm = E (^ cn %- - V , (4) 

1 ^ IS* A fc 



with the coefficients (Z em )^- being pure numbers and D = 4 — 2e standing for the space- 
time dimension. As a result we arrive at the following renormalization group (RG) 
equation for the polarization operator 



' d d \ 

/i 2 ^- + /?(a s )a s — J II = 7 em (a s ) (5) 



or, equivalently, (Lq = In 



u 2 



^-n = 7em (a s )-k)a s i-j n . (6) 
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Here the photon anomalous dimension and the /3(a s )-function are defined in the usual 

way 

7em = /i 2 -^(2 cm ) - eZ cm = + l)(Z em )na l s , (7) 

V?-f^ as = a sP( a s) = - A< +2 - (8) 
a\x j> 

The relation (^j) explicitly demonstrates the main computational advantage of finding first 
the polarization function H(Q 2 ) against a direct calculation of R(s) in the case of massless 
QCDQ. Indeed, in order a™ the derivative ^f^n and, consequently, R(s) depends on the 
very function IT which is multiplied by at least one factor of a s . This means that one 
needs to know II up to order a™ -1 only to unambiguously reconstruct all Q-dependent 
terms in II to order a", provided, of course, the beta function and anomalous dimension 
7e m is known to order a™. On the other hand, the calculation of an anomalous dimension 
or a beta-function is known to be much easier than computing a correlator of the same 
order in the coupling constant (see below). 

On specifying the order of n = 3 we conclude that for computing the O(o? s ) contribu- 
tion to R(s) one needs to know the following ingredients: 



the j3 function to a 2 : it is available through independent calculations of Refs. 



• the polarization function to a 2 including the constant part not depending on logQ 2 : 
the logarithmic contributions were originally obtained in Refs. |24|, |25j while the 
constants were first published in Refs. || albeit not in an explicit form; 

• the photon anomalous dimension 7 em to a^: it is known from works |5|, 0] where it 
was computed with the use of Feynman gauge. 

Thus, in order to check the results of Refs. || || by an independent calculation one should 
compute II in order a 2 s (three- loop diagrams) and 7 em in order a? s (four- loop diagrams). 
The corresponding massless Feynman integrals depend only on one external momentum, 
q in our notation, and will be conventionally referred to as p-integrals in what follows. 



3 Calculation of p-integrals 

It should be stressed that the first, three-loop calculation is definitely much less involved 
than the second, four- loop one. There are three reasons for it. First, the total number 
of diagrams is by order of magnitude less. Second and most important, there exists an 
elaborated algorithm — the method of integration by parts of Ref. |26|, |27j — which allows 



one to analytically evaluate divergent as well as finite parts of any three-loop p-integral. 
Third and also very important fact is that the algorithm has been neatly and reliably| 

1 To our knowledge essentially identical observation was first made in Ref. [Q 

2 The package has been extensively tested in various ways and a probability of the existence of a bug 
in it seems to be extremely small. In fact, the calculation we are describing here may also be considered 
as another highly non-trivial check of the program. 
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implemented in the language of FORM |28| as the package named MINCER in Ref. |29[ . 



An important feature of the algorithm is its ability to evaluate bare dimensionally regu- 
lated diagrams. This allows a convenient two-step evaluation scheme: first calculation of 
bare diagrams followed by a renormalization procedure. The latter eventually reduces to 
a straightforward substitution of the bare coupling constants expressed through the renor- 
malized ones multiplied by proper renormalization constants. Needless to say that such a 
substitution is now routinely done with the help of algebraic manipulation programs like 
FORM. 

The situation with four-loop diagrams is quite different. At present there is simply no 
way to directly compute the divergent part of a bare four-loop p-integral, not speaking 
about its constant part (see in this connection a discussion in Ref. The best what can 

be done is the use of the method of Infrared Rearrangement (IRR). The method is based 
on on an important observation of Ref. |3T| that in the MS-scheme any UV counterterm 



has to be polynomial in momenta and masses. The observation was effectively employed in 
Ref. |^] to simplify considerably the calculation of UV counterterms. The trick essentially 



amounts to an appropriate transformation of the IR structure of Feynman integrals by 
setting zero some external momenta and masses (in some cases after some differentiation 
is performed with respect to the latter). As a result the calculation of UV counterterms 
is much simplified by reducing the problem to evaluating p-integrals. The method of 
IRR was ultimately refined and freed from unessential qualifications with inventing a so- 
called -Reoperation — a generalization of the BPHZ i?-operation to subtract UV and IR 
divergences — in Refs. EBI, . 



The following theorem has been proven in Ref. |34| by the explicit construction of the 
corresponding algorithm: 

Any UV counterterm for any (h+l)-loop Feynman integral can be expressed 
in terms of pole and finite parts of some appropriately constructed (h)-loop 
p-integrals. 

The above theorem coupled with the the integration by parts method solves at least in 
principle the task of analytical evaluation of 7 em to a? s . It should be noted that the 
R* -operation is absolutely essential to prove the Theorem, though in most (but not 
in all) practical cases one could proceed without it. However, such a practice, in fact, 
forces diagram- wise renormalization mode, what, in turns, brings down a heavy penalty 
of manual treatment of hundreds of diagrams. 

Indeed, in genuine four-loop calculations the reduction to three-loop p-integrals is far 
from being trivial and includes a lot of manipulations. Typical steps here are 

a to reduce the initial Feynman integral to logarithmically divergent ones via a proper 
differentiation with respect masses and external momenta; 

b to identify UV and IR divergent subgraphs of the resulting integral; 

c to remove in a recursive way the corresponding UV and IR divergences; 

d to compute resulting p-integrals. 

Among these steps only the calculation of p-integrals can be at present completely per- 
formed by a computer. All others, especially b and c are difficult to computerize. As a 
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result in works [IT], || |6| an inherently time-consuming and error-prone way of manually 
handling every separate diagram (including its full UV renormalization via i?-operation) 
has been employed^. 

In present work we will use the full power of the -Reoperation to simplify the steps 
b and c so far that both UV and IR renormalizations can be done in a global form and, 
consequently, can be simply performed by computer. 



4 IRR in a global form 

We start from the Dayson-Schwinger equation for the correlator ([I]) written in the bare 
form[] 

lC(g, a°) = - J dp^Tr[ la G°(p + q, a°)r° (p, q, a° s )G°(p, a° s )]. (9) 

Here G° and T° are the full quark propagator and the EM current vertex function re- 
spectively; below the integration with respect to the loop momentum p with the weight 
function will not be explicitly displayed. 

The renormalized version of ([]) can be written in two different but eventually equiva- 
lent forms: in terms of renormalization constants or through i?-operation. The first form 
reads 

U aa (q, a.) = Z™q 2 (l — D) — ^Tr[ la G°(p + q, a s )T° a (p, q, a° s )T°(p, a% (10) 

Here Z 2 is the quark wave function renormalization constant; Zy is the renormalization 
constant of the vector current defined as 

where the current inside squared brackets is the renormalized one. The QED Ward 
identity implies the equality Zy = Z 2 , whence the equivalence of (0) and (|T0D follows. 
The second representation is 

U aa (q, a,) = Z cm q 2 (l -D)- R'Tr[ la G°(p + q, a s )T° a (p, q, a s )G°(p, a,)], (11) 

where R' stands for the "incomplete" i?-operation which, when applied to a Feynman 
integral, subtracts only all its UV subdivergences not touching the overall one. 

From the finiteness of the renormalized correlator we have two ways of finding Z em , 

viz. 

Z em = -K e 1 2D{ p _ (^U q Tr[ 1& G\p + g, a°)r° (p, q, a° s )G°(p, a°)] (12) 
+ ^□ g Tr[ 7a G (p + g,a°)r°(p,g,a°)G°(p,a°)]^) (13) 



3 The calculation of Ref. Q did not use the .Reoperation at all while that of Rcf. ||] employed it only 
for a few diagrams. 

4 For simplicity we set the 't Hooft-Veltman unit of mass [i equal to 1 below. 
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and 

Z cm = - K ^{ 2D{ d-1) R> n i Tr ^ G °(P + 9. a ^ (P» 9, a.)G°(p, a,)] | , (14) 
where K e f(e) stands for the singular part of the Laurent expansion of /(e) in e near e = 



and 5Zy = Zy — 1. In Eqs. (|T^) and (14) we have let a Dalambertian with respect to the 
external momenta q act on quadratically divergent diagrams to transform them to the 
logarithmically divergent ones. We also have introduced an auxiliary mass dependence to 
a quark propagator — the one entering into the "left" vertex 7^ — by making the following 
change 

7« -»• 1& = 7« P 2 /(P 2 - ml). (15) 

Note, please, that the auxiliary mass dependence has caused somewhat more compli- 
cated structure of UV renormalizations in the right hand side of Eq. (|T3|). In the case of 
the second equation corresponding modifications are taken into account automatically by 
.R'-operation. 

Now, to IRR. The idea of the method is quite simple: since the renormalization 
constant Z em does not depend on the momentum q dimensionfull one could significantly 
simplify the calculation of Z ern by nullifying the momentum Eqs. ( |i~3|) and Ql4]) . The only 
requirement which must be respected is the absence of any IR singularities in the resulting 
bubble integrals. Unfortunately, a mass introduced to a propagator is not always sufficient 
to suppress all IR divergences. For instance, if q = then there appear completely massless 
tadpoles in the second term on the rhs of Eq. fll3|) . On the other hand, in the diagram- 
wise calculations of Refs. [TT], |5], |J, based on Eq. QT3|), it proves possible to tune the 
position of the massive propagator for any given diagram in such a way to suppress all 
IR singularities. 

Our idea, instead, is to use Eq. (|13|) supplemented by the corresponding IR subtrac- 
tions as prescribed by the -Reoperation formalism [34]]. The problem is facilitated by 



the fact that, as shown in Ref .[35] the IR counterterm constants for a given diagram 
may be determined in terms of some properly chosen combination of UV ones. The only 
remaining task is to write the IR subtractions in a global way. We have done it with the 
following result: 



-K e 



—A— T -O q Tr[ 1& G°(p + q, a s )T° a (p, q, a° s )G°(p, a° s )]\ q = Q 



Several comments are in order regarding this formula. 

First, Eq. ([U|) is, rigorously speaking, applicable as it stands only to the so-called 
non-singlet diagrams, that is to those where both EM currents belong to one and the 
same quark loop. The four singlet diagrams, violating this requirement, appear first in 
order a z s . A full derivation of flTEp with modifications necessary to include singlet case will 
be presented elsewhere. 

Second, by SY\{jp, q, a°) we denote the vertex function of the electromagnetic current 
with the tree contribution removed. The "tilde" atop the index a again means that in 
every diagram the quark propagator entering to the vertex j a is softened at small momenta 
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by means of the auxiliary mass mo according to Eq. fllSD - The bare coupling constant a° 
is to be understood as a s = Z a a s , with Z a being the coupling constant renormalization 
constant. 

Finally, an inspection of (|T^) immediately shows that, in order to find the (n + l)-loop 
correction to Z cm , one needs only to know the renormalization constants Z 2 and Z em to 
order a™ as well as the bare Green functions 

G°(p,a% ^-n(p,q,a s )}\ q = Q , D jr°(p, Q, a°)]\ q = Q , *r°(0,0, a° s ) (17) 

up to (and including) n-loops, that is to order (a°) n . Thus, we have achieved our aim and 
obtained a general formula for Z em in terms of bare p-integrals with explicitly resolved 
UV and IR subtractions. 



5 Results and discussion 



We have computed with the program MINCER the unrenormalized three-loop Green 
functions ([[7]) as well as the quark wave function renormalization constant Z 2 to order 
a%. The calculations have been performed in the general covariant gauge with the gluon 
propagator (g^ — (,^r)/q 2 . We have also taken into account the singlet diagrams as well 
as extra diagrams with some of virtual quarks replaced by colour octet neutral fermions. 
In the minimal supersymmetric standard model such a fermion known as gluino appears 
as the superpartner of the gluon. The total calculational time with an DEC workstation 
exceeds 200 hours for the general gauge; for the Feynman one it is reduced to about 20 
hours. 

Then we have used Eqs. fllED and (]|) to find •y em to order a%. We have used the 
following values for the coefficients of the beta-function in QCD with gluinos 



36 



A) 



1 

16 



jC A - -Tn f - -C A rig 

34 „ 2 , n T 20^ 16 2 

—C A - AC F Tn f - —C A Tn f - —C A n- g 



(18) 



Here C A and Cp are the Casimir operators of the adjoint and quark (defining) repre- 
sentations of the colour group; T is the normalization of the trace of generators of quark 
representation Tr(t a t b ) = T5 ab ; rif is the number of quark flavours; d[R] is the dimension 
of the quark representation of the colour group and rig is the number of neutral colour 
octets which we take either zero or one. 
Our results for 7 em and R(s) read 
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C F Tn 



f 



11 

36 



C F C A n~ q [- 
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7 



23 \ 



32 
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169 11 
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C(3) )+C F C A Tn fi 
V 216 9 v 7 ; V 3888 9 
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<7 
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19 
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C(3) 



(19) 



r 3 



R(s) = d[R]J2Q} 1 + a a (ji)n + a 2 (/x) 
f 1 



s 2 + m^- ( Sl $j) 
s 

S3 + in^l ( 2s2/ g + s^i) + ln 2 ^ ( Sl $ 



(20) 
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128 
90445 



r 2 = C 2 F 



3 

32 



dabcdabc 


■176 


1024 


i 3 




123 


11 




[ 32 





128C(3) 
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+ C F Tn f 



11 



+ C(3) 
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485 11 



55 
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- 3 r- + iHc(3) + |c(5) 



27 72 



C F T 2 n) 



29 19 

-^ + ^C(3)-5C(5) 
54 36 9 sv ; 



+ C 2 F C A n 6 



+ C F C 2 A n- q 



54 36 9 

33767 11 2 251 . 5 . , r , 

1 7T H C(3) H C(5) 

3456 144 36 SV ; 12 sv ; 



+ C F C 2 A n 2 



151 _ J_ 2 _ 19 
216 144 36 U ' 



(21) 



where ci abc = 2Tr({t a t 6 }t c ). 

We observe that neither 7 em no -R(s) depend on the gauge fixing parameter £ as it must 
be. If rig is set to zero then R(s) is in complete agreement with the results of Refs. f| ||]. 
For the standard QCD SU C (3) colour group values C F = 4/3, C A = 3, T = 1/2 and 
^a&c^afcc = 40/3 we get for R(s) with ii 2 = s 



R(s) = 3j2Q}\ l+«s+ 



+ aZ 



87029 121 



365 
~24~ ~ 
1103 
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11 



11 C(3) - - n f + - C(3) n/ - - n § + 2 C(3) 
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7r n/-n 5 r(3)n/-n 5 H n- it n~ n C[o)n- 

18 1 ' !) ; ' 18 9 12 9 •'' 

+ ^ s (j2Qf) fl"^(3)) (22) 



/ 

or, in the numerical form, 



R( s ) = 3 J2 Q}{l + a s + a 2 s (1.98571 -0.115295^- 0.345886™ 

/ 

+ a 3 s (-6.63694 - 1.20013n/ - 0.00518nJ - 2.85053np 



9) 



- 0.03107n/n 5 - 0.04661n?)} 

J2Qf^j 1-2395. (23) 

At last, for the phenomenologically relevant case of rif = 5 we obtain 

R( s ) = — h + a s a 2 (1.409 - 0.346ns) + +a^(-12.805 - 3.006^ - 0.0466n?)l . (24) 
3 L J 

Thus, the 0(ot z s ) gluino contribition to R(s) has the same (negative) sign as in the a 2 
order. We should also add that for a meaningful phenomenological discussion of the 
gluino contribution to R(s) one should also take into account the running of the coupling 
constant in the next-next-to leading order. This requires the knowledge of the gluino 
contribution to the three-loop coefficient fa which, to our knowledge, is not yet available 



in the literature (the purely QCD contribution to (3 2 is known from |37|, |38 . 

To summarize: we have suggested a new convenient way to compute the UV renor- 
malization constant of the correlator of vector quark currents. Our final formula ([16]) 
directly expresses the constant in terms of unrenormalized p-integrals, with all UV and 
IR subtractions being implemented in a global form. The formula is useful in carrying out 



completely automatic calculations. In our previous work [|39] similar formula has been 
obtained for the case of the correlator of scalar currents. 



Using the formula and the FORM version of MINCER j2J| we have computed the 
0(oP s ) correction to R(s) in pQCD including light gluino. In a particular case of the 
standard QCD we have reproduced the result of Refs. [H, El. This gives also an extra 



support for the non-accidental nature of the findings of Ref . H0J . An important feature 



of our calculation was the use of the general covariant gauge, which has allowed us to 
demonstrate for the first time the gauge independence of R(s) at 0(aP s ). 
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